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ABSTRACT 

Working with 108,786 Sloan Digital Sky Survey low redshift galaxies we have 
examined the relation between galaxy mass, metallicity, radius, and star for- 
mation rates primarily in the central portions of galaxies. We subdivided the 
redshift range covered in our sample, 0.07 < z < 0.3, into three narrower red- 
shift bins, and three sets of radial size. We show that for 72% of the galax- 
ies the observed gas metallicities, Z^-, are consistent with (i) a quantitative 
physical relation for star formation through episodic infall of gas of metallic- 
ity Zi = 0.125 X 10^^ ± 1.25 X 10^^; (ii) thorough mixing of infalling and native 
gas before onset of star formation; (iii) a star formation rate (SFR) proportional 
to the 3/2 power of the infalling mass rate, Mf, and (iv) intermittent quiescent 
phases devoid of star formation during which the native gas in a galaxy exhibits 
a characteristic elevated gas metallicity, Zq, dependent on galaxy mass, M*, and 
a characteristic ratio of stellar mass to native mass of gas, Mg. Most if not all our 
star-forming galaxies with < 2.0x lO^^M©, and many with > 2.0x IO^^Mq 
and large radii appear fed by infall. Smaller massive galaxies with high Z^ and 
high star formation rates show more complex behavior. A mean-field-theory toy 
model for the physics of infall accounts for the (SFR) oc Mf^^ relation and per- 
mits us to estimate the mean densities and velocities of clumps of baryonic matter 
traversing the dark matter halos in which the SDSS galaxies may be embedded. 

Subject headings: galaxies: evolution 
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Introduction 



The causes and conditions surrounding star formation appear to be highly complex. 
Despite many attempts to fully understand the different processes that may be involved, 
a globally applicable theory still appears lacking. The past few years, have witnessed 



many attempts to acc ount for star 



or extensions of them (IMadore 



2010 



ormation 



Shi et al. 



aws, s uch as the Schmidt-Kennicutt law. 



20111 ). Until recently, and even now, star 



formation rates have often been treated largely in terms of observa.ble or inferred physical 



conditions within well-isolated galaxies. (ICalzetti &: Kennicutt 



2009 



Krumholz et al. 



20121 ) 



Within the past year or two, however, chemical considerations have begun to invoke 
infall and/or outflow of matter, some of whic h may have observable consequences on 



star formation on ly after considerable delay. (IKobayashi et al. 



Dave et al. 



2007 



Brooks et al. 



2007 



2011ah 



Recently published observational evidence on correlations between star formation 
rates (SFRs) and gas metallicity as a function of galaxy stellar mass is certainly striking 



( iMannucci et al 



20101) • It suggests that star formation may need to be considered in 



entirely new ways that emphasize infall of matter from enveloping halos. 

The Sloan Digital Sky Survey, SDSS, has made available star formation rates and 



(Tremonti et al. 


2004; 


Gallazzi et al. 


2005 




Ellison et al. 


20 


08 




Kewlev & Ellison 


2008 


Mannucci et al. 


2010; 


Pene et al. 


2010 




2011 




Aihara et al. 


I2OII 




Yates et al. 


2OI1I). 



Although these data are largely statistical and may lack the more detailed correlations 
that observations on star formation in nearby galaxies make possible, the advantage of 
working with SDSS data is the ability to establish commonalities in the properties of 
tens of thousands o f galaxies selected from a homogeneously acquired set. In particular. 



Ellison et al 



(120081 ) noted a number of striking correlations between galaxy mass, radius. 
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and metallicity, as well as between mass, star formation rate (SFR) or sp ecific star format ion 



rate ( SSFR), and metallicity. This was complemented by a finding of 



Mannucci et al 



( 120101 ) who produced a three-dimensional Cartesian plot having galaxy mass, star formation 
rate, and metallicity as its orthogonal coordinates. This clearly shows the population of 
star- forming galaxies out to redshifts at least as high as z ~ 0.3 lying on a smoothly curved 
surface. The range of metallicities at each point on this surfa c e is n arrow, revealing the 



su rface to be a thin curv ed sheet. Recent work by 



of 



Mannucci et al. 



Yates et al 



(1201 ll ) analyzes the findings 



(|2010[ ) in greater detail. 



The persuasive correlations these studies have revealed between high star formation 
rates and low metallicity, have made a compelling case for considering infall of pristine, low- 
or zero-met allicity gas as a primary factor governing massive star formation. 

That infall of high- veloc ity clouds is prevalent at least in the Galaxy has recently been 



convincingly documented by 



Lehner fc Howki (1201 ll ) through observations with the Hubble 



Space telescope. They identified the infalling gases through their ionized oxygen, silicon and 
carbon absorption lines seen against stars at high Galactic latitudes, and estimated the rate 
of infalling mass to possibly be of the order of 0.5 Mq yr~^. This is potentially sufficient 
to maintain star formation at an expected rate for the Galaxy. Two questions, however, 
remain unanswered. The observed infall velocities are of order 100 km s~^, which may be 
sufficiently low for the material to simply be gas that is falling back into the galaxy, having 
been previously ejected to some appreciable height. Further, the assumed metallicity of the 
infalling gas, estimated to be Z ~ 0.2Zq, leaves unanswered the question of whether this 
could be the metallicity of infalling gas currently traversing the local halo, or again just 
refiect the metallicity of previously ejected material now falling back into the Galaxy. 



These question s have been given impetus by the work of 



Dave et al 



Dave et al 



( 12011a 



) and 



( 1201 Ibl ). who have emphasized the importance of supernova ejection of 
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high-metallicity material, which then may fall back into a galaxy, depending on whether or 
not the ejection powers escape velocities. 

Evidence for resolving these questions may soon accumulate through observations of 
outflows from star -bursting galaxies, as shown by recent observations on PG1206+459 by 



Tripp et al. 



( l201ll ). While uncertainties about the metallicity of the outflowing gas and 
the extent to which some of this gas may fall back into the observed galaxy still prevail, 
observations of the attempted kind may soon provide informative results. 

Our motivation in pursuing the analyses provided in the present paper has been to 
unambiguou sly quantify the relationship between infall of gas and star formation rates. 



The data of iMannucci et al. 



( I2OIOI ) strongly suggest that the correlation between SFR and 



low metallicity requires the infall of low-metallicity o r pristine ext r agala ctic gas. But other 



Cowie et al.l ( 120 111 ), have not be ruled 



explanations for this correlation, for example those of [; 
out. In these circumstances, it appeared to us that the sheer wealth of data provided by 
the SDSS might serve as a guide, which would permit a more compelling analysis. 

Obtaining an empirical relationship of the type we exhibit in the following sections, 
however, is not enough. Ideally one would like to have a physical model that accounts for 
an observed relation, is consistent with other observational data and, where possible with 
extant theory. For this purpose, we have developed a mean-field-theory toy model that 
indicates how the metallicity and velocity of infalling gas may b e related to SFR and to 
the native gas co nte nt of a galaxy as di s cusse d, for example, by 



Tremonti et al. 



Dalcantonl (120071 ): or 



mm : 



Peeples fc Shankarl (I2OIII ). Such a model should also take into account 



recent calculations on infall of extragalactic gas into galaxy halos, to the extent that 
pristine or low-metallicity gases are most likely to have accumulated in a halo. Moreover, 
the velocity at which this matter traverses the halo and impacts on a galaxy should still 
reflect that history. 
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Although a to y model cannot be a s 



developed, e.g., by 



Brooks et al 



fl2007h : 



comprehensiye a s dedi cate d computer si n iulati ons 



Kobayashi et al 



f l2007h : or 



Klypin et al 



( 1201 ih . it 



provides ready insight on the physical relations between different parameters, and allows 
observers to check for potential confirmation or refutation of hypotheses on the basis of 
readily apparent scaling relations. Within these restrictions, the toy model we develop 
in the present paper appears at least in rough agreement with available data, and makes 
specific assertions that may be negated or verified by future observations. 

A key feature of the present paper is a quantitative demonstration that most, if not 
all star-forming galaxies with < 2.0 x IO^^Mq, and many with > 2.0 x IO^^Mq 
and large radii appear to be fed by infall of pristine or low-metallicity gas. More massive 
galaxies, however, appear to form stars also by other means, possibly through consumption 
of gas injected into interstella r space by outflow from low-mass evolved stars, as suggested 



by 



Leitner fc Kravtsovl (1201 ll ). The evidence we provide is chemical and based on a novel 



quantitative relationship we establish between SFR, metallicity and a galaxy's stellar mass. 

The wealth of galaxies observed in the SDSS permits our toy model to examine not 
only the correlations between such properties as galaxy mass, radius, SFR, and metallicity, 
but, to a limited extent, also the evolution of these quantities with redshift. We restrict 
ourselves to galaxies observed in the redshift range z = 0.07 to 0.3, corresponding to an 
interval of about 2.5 Gyr spanning an epoch that began 3.4 and ended 0.9 Gyr ago. We 
assume standard flat A CDM cosmol ogy, flcM ~ 0.27, Qde ~ 0.73, and a Hubble constant. 



Ho = 71 km si Mpc"^ (IWright 



20061 ). 



We structure our paper as follows: Section 2 details our adopted data processing 
methods. Section 3 characterizes observational results. Section 4 introduces our main 
finding, namely that the relation between metallicity and SFRs in galaxies of fixed mass 
and radius is quantitatively consistent with dilution of native gas in these galaxies by 
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infalling matter. Section 5 discusses our findings in terms of a mean-field-tlieory toy 
model tliat clarifies and quantitatively accounts for the observations. On the basis of the 
SDSS observations, section 6 derives mean values for the velocities {V) and densities (p) 
of infalling matter. In galactic halo models, these appear to be associated with mean 
properties of clumps of gas traversing the halos. In Section 7, we briefly discuss the 
plausibly permitted range of differences between the mass infall rates that trigger star 
formation and the resulting star formation rates. A final section 8 discusses our findings in 
terms of other recently published studies. 



Methods 



Our selection of galaxies comes from the MPA-JHU catalog compiled from SDSS DR7 



available online (http://www.mpa-garching.mpg.de/SDSS/DR7/). 



Mannucci et al 



torn , we 



To permit closer comparison of our data and those of 
have used identical selections in redshift, (between z=0.07 and 0.3), Ha signal to noise 
ratio (>25), an d we similarly exc l uded a ctive galactic nuclei, AGNs, using the criteria 



put forward by 



Kauffmann et al. 



find stellar mass, SFRs, and met allicities. Our ste. 
MPA-JHU catalog as inspired by iKauffmann et al. 



( l2003al ). We have also followed identical procedures to 
l ar mass es are taken directly from the 



from masses based on a Kroupa IMF (IKroupa 



( l2003bl ). with a corre ction to convert 



200l|) to a Chabrier IMF (jChabriei 



Star fo rmation rates are determined using Ha and are based on the work of 



2003 ). 



Kennicutt 



( 1l998al ) with a correction for a Chabrier IMF. 



While Ha is subject to foreground d ust extinction, we cor rect for this using the 



Balmer decrement, along the same lines as 



Cardelli et al 



(Il989[ ). We excluded any sources 



with large extinction corrections corresponding to Ay > 2.5. Metallicities are determined 
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from the strong line diagnostics -R2 3 = ([0//]3727 + 



and [A^//]6584/ifa; as presented in 



Nagao et al 



OUT 



mom and 



4958 + \OIII]5007)/ H13 



Maiolino et al 



fl2008[ ). We 



note that other a l ternat ives for obtaining metallicities have recently been proposed by 



Lara-Lopez et al. 



adopted by 



20101): but for ease of comparisoii we ha ve opted to follow procedures 



Mannucci et al. 



(I2OIOI) and 



Ellison et al 



f!2008h . 



The bulk of our selected sources cover a galactic stellar mass range from ~ 10^ to 
2 X IO^^Mq, with star-formation rates ranging roughly from 0.035 to 6AMq yr~^, and 
oxygen abundances from ~ 4 x 10~^ to 1.4 x 10^^, straddling a Solar System abundance of 
~ 8.5 X 10"'^. Because we concentrate on the additive properties of metallicity, rather than 
metallicity ratios, we express metallicities throughout in terms of actual oxygen abundances, 
as opposed to their logarithmic values more conventionally adopted. 



The redshift range covered both by iMannucci et al.l ( 120101 ) and us is set at a minimum 
value of z = 0.07 to ensure that the [OII]A3727 emission line is well within the useful 
spectral range and that the 3" aperture of the SDSS spectroscopic fiber will sample a 
significant fraction of a galaxy's surface area. At z = 0.07, an aperture of 3" corresponds 
to an angular diameter of 4 kpc implying that we probe only the central 2 kpc regions of a 
galaxy. At z = 0.3, our observations probe galaxi es out to radial distances ~ 6.6 kpc, and 



thus sample the galaxies more fully (IWright 



20061) 



In addition to these selection criteria imposed by iMannucci et al 



( I2OIOI ) and ourselves. 



we required our sample to have an o bserved Petros i an ha lf-light radius in the r-color band 



rsQ. This is required, in the spirit of 



Ellison et al 



( l2008l ). to investigate the role galactic 



radius plays in star formation and metallicity. We have also taken pains to eliminate 
duplicate observations from our SDSS sample. Where an object was observed multiple 
times, we averaged its properties across the multiple entries. 



Both 



Mannucci et al. 



teoioh and 



Ellison et al. 



( 120081 ) have treated galaxies in the 
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redshift range we cover as though they were coeval. To examine whether or not further 
insight could be gained through study of the provenance of their galaxies, we have divided 
the SDSS galaxies into individual tables corresponding to three separate redshift ranges, 
0.07 < 2 < 0.10, 0.10 < 2 < 0.15, and 0.15 < 2 < 0.30; and three Petrosian half-light radii 
rso ranges, small rso < 3.74 kpc, medium 3.74<r5o < 5.01, or large rso >5.01. Each table 
divides galaxies according to logarithmic mass and logarithmic SFR. These tables all are 
available as part of our paper online. 

Although many of the entries in our tables show galaxy populations that may be quite 
sparse in certain galaxy mass/SFR bins, we considered our findings significant only if based 
solely on bins that have > 50 galaxies per bin. 



We are aware that 



Chabrierl (120051 ) has updated his initial mass function and that it is 



in better agreement with current data as indicated by 



Parravano et al. 



order to co mpare our ob s ervat ional data and analytic model with those of 



11). However, in 



Mannucci et al 



feoioh and 



Ellison et al. 



( 12008I ). in section 6 below, we have chosen to adopt an initial mass 
function identical to theirs in order to maintain our respective data sets as homogeneous 
as possible, so that a more direct comparison of our findings and theirs may emerge. This 
procedure admittedly maintains inherent biases that all three data sets may propagate, but 
enables a more meaningful comparison of results. 

While the SDSS provides a way of studying star formation statistic ally in a sample 



exceed ing a hundred thousand galaxies, both our selection and those of 



Mannucci et al. 



teoioh and 



Ellison et al. 



(|2008[ ) restrict themselves to observations of the central portions 



of galaxies, in many of which spiral arms no longer are well defined. Our sample thus is 
not adequate to addressing questions of spiral structure or its effects on star formation, and 
generally obstructs direct comparison to work studying star formation over larger portions 
of galaxies' disks. Nevertheless, given the importance of processes taking place in the 
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nuclear portions of galaxies, our choice of regions has a compelling interest of its own. 



3. The Observed Galaxy Sample 




Fig. 1. — Median metallicity for samples across our range of mass and SFR. This includes 
samples across all redshift cuts and all radius cuts. Here we include only samples with at 
least 50 sources at a given mass and SFR. The small grey points indicate toy model fits to 
our data, as detailed in section 5. Complete numerical tables of these data as well as subsets 
in redshift and galactic radius are available online. 
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In Figure [T] we present our galaxy sample for the entire redshift and radial size range 
divided into logarithmic SFR and stellar mass bins. In each bin where >50 sources exist 
the metallicity is plotted on the vertical axis. 

A full complement of nine sets of tabulations showing galaxy stellar mass, SFR, gas 
metallicity and metallicity scatter, separately for each of the three redshift intervals and 
each of the three ranges of Petrosian radii, is provided online. The tables provide much of 
the statistical information of interest to our paper derivable from the SDSS. 

To detect whether evolutionary effects might have had an impact on previously 
published findings we divided the available data into the three redshift ranges defined 
above. In mass/SFR bins where there were sufficient data in all three redshift epochs our 
tables provide data on the evolution of metallicity with time. There are ~ 40 overlapping 
mass/(SFR) bins in the three red-shift tables, in which each bin is populated by at 
least 50 galaxies. In all, this comprises several thousand galaxies in each redshift range. 
Strikingly, galaxies with identical mass and SFR bins have median metallicities that vary 
across redshift by at most a few percent, somewhat randomly and well within standard 
de viations that g e neral 



of 



Savaglio et al. 



y are of order 15%. This is in satisfactory accord with the findings 
(l2005l ). who observed galaxies in the somewhat higher redshift range, 
0.4 <^ z ^ 1.0 and found no significant evolutionary trends in the galaxy mass-metallicity 
relation. 

Evidently, galaxies with identical masses and star formation rates, in the mass range 
from ~ 7 X 10^ to ~ IO^^Mq, appear to also share other significant properties despite 
the ~ 2.5 Gyr epoch spanning the redshift r ange 0.07 < z < 0.30. In this respect, the 



metallicities exhibited in the single table that 



Mannucci et al. 



(I2OIOI ) provided appear to be 



robust even though they included data on galaxies over the entire z = 0.07 — 0.30 redshift 
range and thus galaxies observed at epochs that may differ by as much as 2.5 Gyr within 



Table 1: Main Sequence Galaxies: Peak Populations in the Mass/SFR Tables 



Redshift Coverage of Table 0.3> z> 0.15 0.15 > z> 0.10 0.10 > z> 0.07 

Metallicity Range Spanning all Masses 

(N(0)/N(H) xlO^) 0.61 - 1.37 0.51 - 1.39 0.473 - 1.34 

Galaxy Mass, M*, in Peak Population Bin 

in units of IO^^Mq 5.62 2.82 1.41 

Metallicity in Peak Population Bin 

(N(0)/N(H) xlO^) 1.26 1.23 1.17 

SFR in Peak Population Bin 

(in units of M© yr-i) 2.239 1.122 0.398 

r = M*/SFR at Peak Population 

(in units of 10^° yr) 2.5 2.5 3.5 



individual bins. 



We find that even at the highest SFRs, the galaxies in our sample could not have 
gained sufficient mass during the ~ 2.5 Gyr between epochs z = 0.3 and 0.07 to shift 
significant numbers of galaxies from one mass column into a neighboring column. As table 
1 indicates, SFRs diminish toward lower galaxy mass but in such a way as to apparently 
enable all of the galaxies to accumulate mass at a rate proportional to their own masses. In 
effect, the specific star formation rate (SSFR), defined as the star formati o n rat e per galaxy 



Peng et al. 



fl2010l ) have also 



mass, remains constant for this peak population component, 
noted that the deviation from constancy of the SSFR is quite small across all masses in the 
SDSS. Here we are finding a related result, namely that their relation holds even though a 
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significant number of liiglier-mass galaxies are observed at an epocli 2.5 Gyr earlier tlian 
virtually all low-mass galaxies. The time galaxies at all masses would require to double 
their mass at their observed SFRs is practically constant, is the reciprocal of SSFR, and 
assumes values of ~ 2.5 — 3.5 x 10 ^° yr, roughly 2 times the age of the Universe at each 



of the observed epochs ( 1 Wright 



20061). 



These data should, however, not be interpreted as indicating low, steady infall rates 
that maintain a galaxy in a relatively stable state of mass accretion. We envisage an 
individual galaxy of given mass periodically moving down or up within its mass column in 
one of our tables, as infall from its surrounding halo increases or decreases both the infall 
and the star formation rate within the galaxy. As we begin to show, starting with the next 
section, mass infall rates into galaxies vary across a range spanning nearly two orders of 
magnitude, need not be constant in time, and appear to be consistent with episodic star 
formation in the SDSS sample we have studied. 



4. Infall and Metallicity Dilution of Native Gas 

We now examine the extent to which the amount of mass MjTj, which has fallen into a 
galaxy before onset of star formation, dilutes the initial metallicity Zq of the mass of native 
gas, Mg, present before the onset of star formation. 

If we assume that the infalling matter has metallicity Zj, and mixes thoroughly with 
the galaxy's native gas, we expect that the observed metallicity should be given by the 
ratio of the aggregated mass in metals after infall, divided by the total aggregated gas. 

^ _ (M,rZ, + MgZo) 

From equation ([1]) we can see that in galaxies which have an amount of infalling gas equal to 
the native gas mass, MjT = M^,, the observed metallicity will be half the native metallicity, 



= |Zo, if the infalling gas has zero metalhcity, = 0. 

To test equation ([T]) broadly, we adopt a trial function expressing mass infall rates, Mi, 
for which we lack observations, in terms of SFRs for which we do have observations. 

Mi = a{SFRy/\ equivalent to (SFR) = eMf^^ . (2) 

Here a and e = a^^^"^ are constants of proportionality having some as-yet-unknown 
value but units, respectively of Mq^ yr"-*^/^ and Mq^^'^ yr-*^/^. 

We will justify the Mi - SFR relations of equation ([2]) in section 5. The second of 
these relationships is the more physical. The factor e is a measure of star formation 
efficiency. But, for now, we will make use of the relationship based on a, because (SFR) 
is an observationally determined quantity, whereas Mi is derived. Using this relation the 
observed gas metallicity takes the form 
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Table 2: Star Formation Rates and Observed and Computed Metallicities for low redshift, 
small half light radius galaxies with M* =O.5xlO^'^M0 in Figure [2j 



(SFR) 


(SFR)^/'^ 


10^ X Zx Observed 


10^ X Modeled 


Mo/yr 




N(0)/N(H) 


(Eq[3]) 


0.10 


0.22 


1.00 


0.99 


0.14 


0.27 


0.97 


0.97 


0.20 


0.34 


0.94 


0.95 


0.29 


0.43 


0.90 


0.92 


0.40 


0.55 


0.88 


0.89 


0.57 


0.69 


0.86 


0.86 


0.81 


0.87 


0.82 


0.82 


1.14 


1.09 


0.79 


0.77 


1.61 


1.37 


0.72 


0.72 



As an example, we may use the entries for low redshift (0.07 < z < 0.10,) medium 
radius (3.74 < r^Q < 5.01 kpc,) galaxies of mass 0.5 x IO^^Mq from our tables. Adopting 
the observed metallicities in units of oxygen abundances (10^xN(O)/N(H)), at the available 
SFRs (Mq yr^-*^), we can calculate best fit values of Zq and {Tia)/Mg. For these data we 
find a good fit with Zo=1.06xlO-3 N(0)/N(H) and {Tia)/Mg = 0.414 (Mo/year) "^/s. Table 
|5] shows the results of this sample calculation and Figure [2] shows equation fitted to 
the sample data. Regrettably, we cannot isolate Mg from axj with the data given, and can 
derive solely their ratio {aTi)/Mg. 

Online we provide our metallicity fits like those of Figure [2] for all stellar masses, radii, 
and redshifts for which we have adequate statistics — by which we mean > 50 galaxies per 
mass/SFR bin and 5 of such bins per galaxy mass column. For each sample of galaxies 
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0.0 0.5 1.0 1.5 2.0 

SFR (M3„„/yeor) 

Fig. 2. — Observed gas metallicities at various SFRs for M* = 0.5 x lO^^M© galax- 
ies with low redshifts and small half light radii. The fitted equation ([3]) with = 
IMAx 10-^ N{0)/N{H), {Tia)/Mg = 0.414 (M0/year)-2/3, and Zi = 0.125 x lO^^ is plotted 
as a solid line. The fit is acceptable with a value of 0.35. Dashed and dotted lines indicate 
fits based on our lower and rough upper limits on Zi, and 0.25 x 10"'^ respectively. Similar 
plots for galaxies with different stellar masses, radii and redshifts are available online. 

we performed a analysis, testing equation ([2]) with a range of 0.1 < x 10^ < 1.5 
N(0)/N(H) and < < 5 (M0/year)-2/3 and selecting the lowest reduced fit- 
Observational uncertainties in the emission lines used in our strong line diagnostics are 
propagated through to find the uncertainty in metallicity for each mass bin. Horizontal 
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error bars represent the uncertainty in SFR due to the SFR bin width. To make use of 
the SFR uncertainties in our analysis we found how much the model metallicity varies 
over the range of each SFR bin and took half this variance to add in quadrature with the 
metallicity uncertainty. This conversion of uncertainty in SFR to uncertainty in metallicity 
is performed for each SFR bin in each galaxy sample with every trial of model parameters. 

While we expect and fp- to vary between galaxies of different mass, in our naive 



picture the infalling gas should not anticipate what kin d of galaxy i t is fa 



should not vary systematically between mass bins. As iDave et al.l ( j2011a( ) and 



ling i nto and thus 



Dave et al. 



(l2011bl ) have pointed out, however, feedback in the form of supernova ejection may pollute 
extragalactic gas, raising the possibility that infalling gas might already have significant 
metallicity. For that reason, in our analysis we allow the infalling gas metallicity, 
Zj, to be nonzero but fix it at the same value for each sample. Repeating the whole 
fitting procedure many times with different fixed values of Zi we find that a value of 
Zi = 0.125 X 10^'^ minimizes the total of combined fits. A range of low metallicity 
solutions, spanning < Zi < 0.20 x 10^'^ provide nearly equally good fits, while for 
Zi > 0.30 X 10~^ the fits become clearly worse and the total value sharply increases. We 
therefore take 0.125 x 10^^ as our infalling gas metallicity with the caveat that our results 
are consistent with any very low metallicity infalling gas Zi ~ 0.125 x 10^'^ ± 0.125 x 10^'^. 

Our metallicity fits show that 72% of the SDSS galaxy samples obey the metallicity / 
star- formation relations and are well fit by equation ([3]). Given our small number of degrees 
of freedom (which ranges between 3 and 13) for each of our metal dilution samples, we used 
a reduced value of 2 for our cutoff. 

For the curves that provide good fits, we averaged both the Zq and Mg values obtained 
for each group of galaxies with a given stellar mass M^,. For the mass sequence from 
= 0.36 to 11.2 X 10^*^ Mq for which we have sufficient data, we find a monotonic increase 
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in metallicity going from Zq = 0.90 to 1.37 x 10^ N(0)/N(H). We normalize the other fitting 
parameter by multiplying it by stellar mass, ariM^,/ Mg. This parameter stays relatively 
constant over the entire galaxy mass range, hovering around an average value of (2.3 ± 0.7) 
xlO^ (yr^ Mq)^/^. On the assumption that a and Ti are constants across our sample bins, 
this parameter is indicative of the stellar to gas mass ratio. If the gas fraction in a typical 
galaxy were Mg/M,, ~ x, the corresponding fitted parameter would be arj ~ 2.3 x lO^x (yr^ 
Me)V3. 

Sets of galaxies not obeying equations ([T]) and ([3]), invariably exhibit high metallicities 
at high SFRs and masses M^, > 2.0 x IO^^Mq, suggesting that for these galaxies star 
formation independent of pristine or low-metallicity mass infall plays a significant role. 

The good quality of the fits we observe convince us that the dilution we have modeled 
actually takes place. Particularly interesting then is that the anticipated dilution must take 
place before star formation is initiated, so that the observed HII regions, which enable us to 
determine SFRs and metallicity values, reflect the conditions that triggered the formation 
of the stars. 



Two points m a y sti. 



1 be noted. 



Tremonti et al. 



t00i \. balcantonl (120071 ). and 



Peeples fc Shankarl ( 1201 ll ) have evaluated a number of trends of gas mass and metallicity 



in galaxies with differ ent stellar masses. However, as the detailed atlas of galaxies compiled 



by 



Leroy et al 



(|2008[ ) makes clear, global ratios of gas-to-stellar mass vary enormously, 
depending on whether one is dealing with a narrowly defined central region of a galaxy, as 
defined by the SDSS s pectroscopic apertu re used throughout the present paper, or entire 



galaxies, as studied by 



Leroy et al. 



(120081 ). Extended galaxies generally show the ratio of 
gas to stellar mass in a galaxy's central regions to be of the order of a few percent, even 
when that ratio, averaged over the entire galaxy is one or two magnitudes higher. But even 
in these central regions, the ratio M^/M* can vary greatly. 
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For this reason, legitimate questions may be raised about our assumption that some 
well-defined metallicity and fractional amount of native gas can be defined for galaxies in a 
given stellar mass range, and that both parameters may be expected to vary systematically 
and continuously over galaxies having different stellar masses. In particular, one may ask 
what physical basis might account for such a relationship. 

The response to these questions is that, in the present investigation, both the gas 
fractions and the metallicities involved are statistically determined. For each of the galaxy 
masses and SFRs considered we deal with groups no smaller than 50 galaxies. For galaxies in 
each stellar mass range, both the derived native gas metallicity and an assumed gas fraction 
derived from all the distinct SFR ranges are averaged over hundreds if not thousands of 
galaxies. To the extent that the SDSS galaxies fall into different classes, governed by 
distinct physical conditions, there may be no single physical reason for assuming that such 
average values actually exist, other than that the occurrence of each class of galaxy in the 
general distribution of galaxies follows some well-defined evolutionary history, and that 
the physical parameters governing these galaxies' gaseous components fall into a relatively 
restricted range. 

Table [3] provides additional insight offered by an analysis of galaxies of different radii. 
In this table we show the fitted native gas metallicity values, Zq, at all masses and radii 
where adequate statistics were available in a mass range spanning 0.25 - 4 xlO^° Mq. Note 
that satisfactory fits using equation (|3]) were found for all three radii, at low and medium 
redshifts at galaxy masses 0.7 < MJ{10^^ Mq) < 1.4 and 1.0 < M,/(10^° Mq) < 2.8, 
respectively. Below this range there were not enough data for a significant fit; above this 
mass range, at the radii indicated, the data were not well fitted by equation ([3]) and are 
indicated by blanks in Table 4. All the galaxy bins displayed show a trend of increasing 
metallicity with mass and generally decreasing metallicity with radial size. 
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Table 3: Native Gas Metallicities (Zq x 10'^) as a Function of Stellar Mass and Galaxy Radius. 
(For details see text.) 





0.10 > 


z > 0.07 


0.15> z > 0.10 


0.3> z> 0.15 




r5o: Small 




T , f1 T* CTfi 


^Tn 51 1 1 

Oliidil 


iV±CJ.. 


T TCffi 






0.36 


0.90 
















0.50 


1.06 


l.UU 




"I '^f\ 

I.OD 










0.71 


1.18 


1.13 


1.06 


1.32 










1.00 


1.25 


1.17 


1.18 


1.30 


1.20 


1.14 






1.41 


1.30 


1.24 


1.24 


1.39 


1.28 


1.22 






2.00 








1.39 


1.33 


1.24 






2.82 








1.37 


1.33 


1.30 




1.30 


3.98 






1.31 


1.35 


1.34 




1.41 


1.29 


5.62 






1.27 


1.34 






1.31 


1.31 


7.94 










1.37 


1.30 


1.36 


1.31 


11.2 
















1.37 



The significance of the parameter a remains to be determined. This becomes apparent 
if equation ([2]) is recast in the equivalent form (SFR) = eMf^'^, where e is a constant whose 
units are (M0/?/r) ~^/^and a = e~'^^^. Seen in this form, e is readily interpreted as a factor 
determining star formation efficiency and inversely proportional to the cube root of the 
mass infall rate. In section 5, immediately below, we will show how this efficiency relates to 
pressure regulated star formation in triggered molecular cloud collapse. 

Before leaving this section, we should still address the question of mass balance. Gas 
falling into a galaxy should either add to the galaxy's gaseous or stellar mass, or add mass to 
a central black hole, unless an equal amount of mass is ejected from the galaxy explosively 
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or by means of a variety of potential winds. Equation ([2]) thus should not be interpreted 
as meaning that the duration of infall is identical to the duration of star formation. As 
discussed in the next several sections, we deduce that the two processes are sequential, and 
that their respective durations are generally not identical. 



5. A Mean-Field-Theory Toy Model 

To account the findings reported thus far within the larger framework of astrophysics, 
we next adopt a toy model based on a mean-field-theoretic approach. 

Consider a distant, highly tenuous, cold clump of gas of mean mass density (p) directly 
approaching a galaxy of mass M. Let the relative velocity of the clump at large separations 
be V . This clump may have been previously captured into a massive halo surrounding the 
galaxy. We may follow the motion of a volume element of the clump displaced by an impact 
parameter s from the point on the clump that will eventually impact on the galaxy's center. 
At large distances from the galaxy, the angular momentum of a unit-mass increment of gas 
about the galaxy's center is sV . Because gas pressures within a cold tenuous clump may be 
largely neglected, angular momentum conservation tells us that the velocity v of the mass 
increment impacting on the galaxy at radius r is given by 



Vs = vr , (4) 
and the infall energy of the increment is 



V MG 

Here G is the gravitational constant and, as usual, the two equations deal solely with 
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the dominant azimuthal velocity components while neglecting minor helical velocity twists. 

Within some radial distance from the galaxy's center, infalling matter will be 
captured through collision with native gas gravitationally anchored to the galaxy. Beyond 
Tm the galaxy's gaseous component may be too tenuous and/or M may be too low, and 
capture will cease. Solving the two equations for s and V we obtain a capture cross section 
for matter falling into the galaxy within radius r^. 



TTS 



■ 2 2rmMG 
r H 



(6) 



Within a radial distance from a galaxy's center, the mean infall rate depends on the 
mean intergalactic gas density (p) as 



(M,) = 7r(p)sV = 7r(p) 



■ 2 , 2r„MG' 
r H 



V . 



(7) 



Every indication given by the low metallicity of SDSS star-forming galaxies is that star 
formation at least in these galaxies is initiated by infall of low-metallicity gas. Triggered 
star formation is ge nerally thought to com mence when the pressure on a molecular cloud 



is sufficiently high (IKrumholz et al. 



20091 ). The pressure P exerted on the galaxy due to 
infalling gas will vary as a function of radial distance r, and will also depend on infall angle, 
i.e., whether impact of the infalling gas is at a glancing angle or perpendicular to the gas 
layer it impacts. But because our data provide only a single star formation rate per galaxy, 
the parameter of immediate relevance may be some representative pressure (P) suitably 
averaged over radial distances r and angles of incidence 9. Taking v now to be the infall 
velocity at radial distance r, and v{f{6)g{r)) a representative perpendicular-impact-velocity 
component averaged over all radial distances r and directions 6, we can write the pressure 
for infall onto an area A within radial distance from a galaxy's center as 
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(P) 



Miv{f{e)g{r)) 
A 



(P) 



7rr, 



2 r 



A 



1 + 



2MG 



3/2 



(8) 



where we have made use of equations 



(ISl) and (ED to ehminate v. 



An important feature of equation is that (P) is solely a function of the infall rate 
Mj, galaxy mass M, and radius r^, as well as two universal parameters characterizing 
the infalling matter, a representative velocity for clumps traversing the halo, V ^ and a 
representative clump mass density, (p). Although we have not yet derived estimates for V 
and (p) on the basis of available data, we will do so in section 6, below. 

Massive star formation is known to occur principally, if not solely, in giant molecular 
clouds. Two recent findings regarding these clouds now enter consideration: 



Krumholz et al.l (|2012[ ) have shown convincingly that star formation rates are directly 
proportional p^^^, where Pct is the molecular cloud density, in the Galaxy, local disk galaxies 
and star bursts, and high redshift disks and starbursts. The universality of this relationship 
is documented by their Figure 3, which establishes a single observationally derived ratio 
of SFR-per- unit-area S^, to areal-density-divided-by- free-fall-time, S/t/j. Within factors 
of order ±3, this ratio remains constant at ~ lO^^S/tjj over a SFR density range 
of six orders of magnitude. Here, S is the areal density, t// = A/37r/32G'p^^, and G is 
the gravitational constant. Two relevant factors to take away from this is that, both on 
dynamic and observational grounds, (SFR) oc p^^^^, and that the observed star forming 
efficiency is of order 1%. 



Juvela fc Ysardl (1201 ll ) show theoretical 



A second consideration is that calculations by 
considerations to imply molecular cloud core temperatures to be remarkably constant. 
Equilibrium temperatures, at which cosmic ray heating is balanced by spectral line cooling 
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and cooling by dust to which the gas is coupled, are expected to cluster around T ~ 10 
K for clouds in the density range n{H2) = 3 x 10^ to 10^ c ni~^. These coii s i derat ions are 



Solomon et al. 



(119871 ). Gas 



in general accord with early temperature measurements by 
clouds evidently undergo rapid cooling down to temperatures of order lOK, below which the 
cooling becomes less efficient. The temperature constancy across such a wide density range 
suggests that, over sufficiently long periods to permit heat to be radiated away, molecular 
clouds thus may be isothermally compressed through application of an external pressure. 

We now examine whether gas falling into a galaxy can produce the required pressure 
to significantly increase the density of a molecular cloud. Using a virial coefficient defined 

= 5af^^R/MgG, where R is the molecular cloud radius, Mg is its mass, and 
atnt is the one-d i mensi onal thermal plus turbulent velocity dispersion within the cloud. 
Krumholz et al.l (l2009l ) give the mean pressure in a cloud as Pee = 0.7 x lO^aMrS^/ ks dyn 
cm~^. Here S^^ is the areal density of a cloud in units of ~ 85Mq pc~^ ~ 2 x 10~^ g 
cm^, ks i s the Boltzmann con stant, and a^r ^ 2 independent of ga laxy environment, as 



found by 



Bolatto et al. 



(120081 ) and revised by 



Krumholz et al. 



( 120091 ). A value of 



implies that the cloud is in self-gravitating virial equilibrium, so that the internal turbulent 
plus thermal pressures approximately balance gravity. When amr ^ 1, a cloud undergoes 
free-fall collapse unless magnetic fields support it. 

To exceed a cloud virial coefficient a^ir ~ 1 an external pressure needs to typically 
exceed P ~ 10^^^ dyn cm~^. We can use equation ([8]) to determine the conditions under 
which pressures of this magnitude may be expected. Setting the factor {f{6)g{r-)) ~ 1, we 
then obtain 



P = 10 



-11 



3/2 



O.SMo/yr 



V 



170 km/s 



^13 kpc^\ /3.3 kpc\ / 5 X lO^^s g/c^i^^ 



A 



dyne / cm^ 



(9) 



where we have anticipated physical parameters we will justify in section 6 below. The point 
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to recognize, however, is that an infall rate of ~ 0.5Mq per year converging on a galactic 
area smaller than ~ 13 kpc^ will generate a sufficiently high pressure to compress typical 
galactic clouds and thus trigger star formation. 



The trig 



g ering 



Harwit et al. 



process no doubt is complex. Early calculations of supersonic infall by 



(119871 ) showed that, despite a succession of layers of highly or partially ionized 
gases, and neutral components separating the tenuous, high-velocity impacting gas from 
a cold, relatively dense neutral layer in a galaxy's plane, the pressure P remains constant 
throughout such a layered stack. Thus the pressure exerted on a galaxy's impact area, 
if large compared to the scale height of the gaseous component, is also the compressive 
pressure on a cold cloud in the galaxy's plane. Recent studies devoted to infall into 
massive cluster haloes rather than sing le galaxies are based on similar hydrodynamic 



co nsiderations ( 



by 



Birnboim &: Dekel 



Faucher-Giguere et al. 



20111 1. Roughly similar approaches have been undertaken 
( 1201 ll ). While highly interesting, these studies still await 
observational confirmation, which may take some years to obtain. The connection to the 
infall of clumps into haloes, and their ultimate impact onto the individual galaxies observed 
in the SDSS survey, will be discussed below, in sections 7 and 8. 

Supersonic infall creates turbulence that intersperses infalling gas with cold neutral 
components and leads both to compression and eventual diffusive mixing. The rapid cooling 
of molecular gases, implied by the uniformly low 10 K temperatures pointed out earlier, 
ultimately also dissipates turbulent structures, leaving the cold gas compressed, at the same 
temperature T ~ 10 K it had before infall, and with a sufficiently low virial coefficient to 
collapse. The compression ratio \E' in this adiabatic sequence is ^ 4; its value will depend 
in part on the ratio of heat capacities 7 = Cp/ c„ in the various turbulent gaseous strata, 
and in part on the rate of conversion of atomic into molecular hydrogen. Even a value of 
\I' = 4, however, will diminish the free-fall collapse time by a factor of v^, and thereby 
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greatly accelerate, even if it does not initiate, collapse and star formation. An acceleration 
by a factor of 3 in a pre-stellar free-fall collapse, which otherwise might require ^ 10^ yr is 
significant. Galactic rotational periods are of order 10*^ yr, so that impact of extragalactic 
matter on a particular area of a galaxy, will generally produce pressure on this area for 
no more than ~ 10^ yr. Over longer periods, the galaxy's rotation gradually moves an 
impacted area out of the infalling stream that subsequently may fall on neighboring areas. 

In summary, clumps of intergalactic matter in haloes, in the density range and with 
a mass-infall rate indicated by equation appear able to both shock-compress and 
turbulently mix with native cold gases in a galaxy, in times that suffice for the gas to cool 
to ~ 10 K, dissipate the turbulence, and collapse in free-fall over times of order 10'' yr. 



Elmegreenl ( 1l993l ). which lays the theoretical 



Relevant to the discussion also is work by 
foundation linking the pressure of the ISM, and the ambient UV radiation field to the 
resulting fraction of hydrogen in its molecular form, predicting a nearly linear relationship, 
Kmoi ocP^'^ where Kmoi represents the fraction of hydrogen in molecular form. Recent 



work by other authors have refined this prediction based on empirica. 



simulations, lowering the power law index to ~0.9 (jWong fc Blitz 



2006 



Robertson fc Kravtsov 



2008 



Leroy et al. 



20081). 



2002 



observations and 



Blitz fc Rosolowsky 



None of this is meant to deny that supernovae trigger star formation as well. The 
observational evidence for their propagation of star formation is clear. But the strong 
correlation of low metallicity with star formation in SDSS galaxies suggests that star 
formation is initiated by infall; once initiated, star formation is then propagated further 
by supernova action. Hydrodynamic instabilities have a role as well, but the additional 
pressures exerted by infall and/or explosive ejecta can compress and accelerate the collapse 
of unstable neighboring regions to give rise to observed O and B star associations. In 
addition, although it is by now well established that regions of star formation line spiral 
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arms, the limited spectroscopic aperture of the present study focuses attention on the 
central portions of galaxies, where spiral arms become less distinct. Our data thus lack 
observations that could meaningfully add to a discussion of star formation in spiral arms. 



Finally, we note that lEUison et al.l ( l2008l ) originally demonstrated a correlation of 
higher metallicity with lower radii in galaxies of identical mass, and a similar correlation 
of higher metallicity with lower star formation rates. Equation ([7]) provides a physical 
explanation of why infall rates should diminish with diminishing radius; and equation ([1]) 
shows why lower infall rates diminish dilution of the intrinsic metallicities of native gas 
clouds and thus lead to the observed higher metallicities. Equation ([2]) further shows that 
diminishing infall rates result in diminishing star formation. 



It may be useful, at this point, to still mention that the relation between (SFR) oc Ml 
in equation ([2]) has a superficia. 



V3/2 



T^^sFR) OC T^l'^^ dKennicutt 



resemblance to the Schmidt-Kennicutt (SK) law, 
1998bl ). Roughly equally good fits to the metallicity dilution 
could be found with (SFR) oc M^-*. The SK law, however, deals with the areal density of 
the native gas in a galaxy, whereas lilf^'^ refers to the pressure infalling gas may exert. 



We next turn to a justification of the physical parameters given in equation (Q. 



6. Representative Values of s and V 

The wealth of available SDSS data permits us to derive representative values of s, V 
and (p) by means of relations IQ, ([7]), and ([2]). We obtain these values for a sampling of 
galaxies on the galaxy main sequence, i.e., for the most populous and typical galaxy sample 
available in our data. We define our observed main sequence as the straight line in log-log 
space: Log(SFR)=1.12xLog(M*)-11.74. Along this straight line we have picked five sets of 
values of galaxy mass and SFR more or less equidistant along the most populated portion 
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Table 4: SFRs, Masses, and Radii for Typical Main Sequence Galaxies 



Galaxy 


(SFR) 


/ j-N -1—1 -I— » \ o /q 

(SFR)2/^ 


M* 


^^50 


# 


(p)/a at V=170 km/s 


Set 


Mo/yr 


(Mo/yr)2/3 


Mo 


kpc 


sources 


(g cm"3)x(s g^-^y^^ 


1 


0.209 


0.352 


0.71 X 10^° 


4.2 


62 


1.0x10-36 


2 


0.417 


0.558 


1.3 X 10^° 


4.1 


141 


1.3x10-36 


3 


0.813 


0.871 


2.4 X 10^° 


4.4 


167 


1.3x10-36 


4 


1.51 


1.32 


4.0 X 10^° 


4.9 


127 


1.3x10-36 


5 


3.16 


2.15 


7.9 X 10^° 


5.6 


64 


1.0x10-36 



of our main sequence. At each SFR/M* position we extract a set of sources within a small 
Log(SFR) / Log(M*) range of ±0.02 and take the median Petrosian half-light radius to be 
representative of the entire set. The five sets are summarized in Table HI 

We do not know precisely where the infalling material is stopped through interaction 
with gas in the galaxy, or exactly what fraction of the mass within that radius is stellar, as 
contrasted to dark matter. We assume in the equations below that infalling mass is stopped 
at around two Petrosian radii, and that the enclosed mass within that radius is composed 
largely of stellar and da rk matter w i th a t otal mass of ~ 2xM*. This is roughly consistent 



with data published by 



Leier et al. 



(120 111 ), who plot ted the dist r ibutio n of stellar mass 



Leier et al. 



( 1201 ih Figure 5). Their 



and total gravitational mass in lensing galaxies (see 
total enclosed mass, which they term lensing mass, plotted as a function of effective, i.e., 
half-light radius, varies with galaxy stellar mass, but for the range of masses we consider, 
their ratio of total mass to stellar mass, enclosed within two half-light radii, is roughly 2 to 
3M^, with sizable error bars. 
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We next solve for p and a using equations ((6]), ([7]), and ([2]), finding 

1 SFR^/^ ( 2MGY^ , , 

Based on the assumptions that V has some common characteristic value and (p) is 
also more or less constant when averaged over infall rates on main sequence galaxies, we 
consider a range of V values and examine the resulting (p) / a values for each of our five sets 
of sources. Due to scatter in the data, some difference in the resulting (p) / a values for our 
five sets is inevitable. In figure |31 however, we show that there is a distinct minimum in 
this scatter at a V value of 170 km/s. Specifically, this is where the standard deviation of 
{p)/a among the five source sets normalized by the mean value, reaches a minimum. The 
corresponding mean (p)/« value is 1.2x10"^^ (g cm~'^)(s g^^)^/^ . 

That the infall velocity at large distances, V is as high as ~ 170 km s~^ may be an 
indication that gas flowing into a halo containing a group of galaxies acquires this high 
velocity on falling into the halo, and then impinges on the galaxies without significant 
further acceleration. For a halo with radius ~ 300 kpc and mass ~ lO^^M©, infall velocities 
of this order would be expected; a more precise value would require knowing the radial 
distribution of the dark matter within the halo. Velocities of this order may thus be 
compared to those obtained from simulations of infall into haloes such as those recently 



calculated by 



Klypin et al. 



mm . 



Our estimate of an approach velocity of 170 km s~^ is in rough agreement with 
estimates of V based on the dispersion of line of sight velocities of galaxies in groups and 
poor clusters, with a range of radial (line-of-sight) veloci ty dispersion of 10 0-500 km s~^ 



and a median radial velocity dispersion of ~ 250 km s ^ (lAUen fc Cox 



19991 ) 



Adopting a representative mutual approach velocity ~ 170 km s for present 
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Fig. 3. — Scatter of p/a function of V 

purposes and setting a{SFRy/^ equal to the mass infall rate Mj,we solve for (p) with the 
help of equations (j7]) and ([2]) , by writing 

Mi a{SFRf'^ 



ip) 



(11) 



With the parameters given for the five galaxies in Table IH we obtain a mean mass 
density (p) = {^^^^-^)4.8 x IQ-^^ g cm'^. 



In section 7, below, we seek a value of a consistent with available observations in order 
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to obtain an actual estimate of mean clump density (p). 

7. Difference Between Infall and Star Formation Rates, a. 

In section 4, we found that aTiM^/ Mg has a median value of ~ 2.3 x 10^ (yr^ M©)"*^/^. 
For a typical galaxy with M^/M* = 5%, this yields ari ~ 1.15 x 10^ (yr^ M©)-*^/^. We now 
return to our example of Figure [21 where we dealt with galaxies of mass 0.5 x lO^^M©. It 
will be useful to recast equation ([T]) in the form 

M,r, = aSFR^'^n = ^° ~ ^" M,. (12) 

Looking at the metallicities of the most vigorously star forming samples in Figure [21 
we find that the infalling mass of gas corresponds roughly to 57% of Mg. For MgjM^ = 0.05 
we then obtain an infall time of ~ 1.4 x 10^MQ/(a(SFRf/^) ~ 1 x lOV ( ^j^^ — - ) 
yr. A value of a = 10, would reduce the required infall time to ~ 10^ yr, comparable to 
the time that gas falling onto a limited portion of a galaxy would impinge there before the 
galaxy's rotation moved neighboring regions in the galaxy into the infalling stream. 

An alternative estimate of a is obtained from parameters derived in section 5 that 
would produce a pressure sufficiently high to trigger or accelerate free-fall collapse. 
For P ~ (Miv/A) ~ a{SFRy^^v / A and an infall velocity, as given by equation (jS]), 
V = {s/r)V ~ V3V ~ 294 km s-\ the infall pressure P{IMq yr"^ kpc"^) ~ 1.3 x 10"^° 
dyn cm~^, where we have taken \^ = 170 km s^^. Exceeding the internal molecular cloud 
pressure then would require a > 1. 

Another estimate of a may be suggestive: Spectra have by now been obtained for 
approximately 6 x 10^ SDSS galaxies in our redshift range. Roughly 17% of these have 
indications of star formation and are included in our tables. If 17% of the galaxies have 
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measurable star formation rates, 17% of a typical halo volume should be filled with clumps 



of density (p) = 4.8 x 10 ( ^^^/g" ) g cm ^, and the gaseous mass density of the haloes 



should be 8.2 x 10 ( — tt^ — tt: I S cm ^. For the above-cited halo of radius 300 kpc, the 



total clump mass would then be ~ 2.7 x 10^^ ( j^^i/^'^^.^/^ ) g, or 13.7 x W^{jj^^)M^ 



LQ. 



This suggests that a should not greatly exceed a ~ 1 to 3; otherwise the fraction of clump 
mass in the halo would begin to approach that of the total halo mass estimated above as 

In summary, it is possible to defend a values in the range of 1-10 (Mq yr~^)^/^, 
representing the difference between values of {SFRy^^ and Mj. This also appears roughly 
compatible with respective durations of Ha emission from HII regions around massive stars 
formed in the collapse of a giant molecular cloud and the duration of gas infall onto a 
particular spot of a rotating galaxy. 

Ultimately, however, the value of a will have to be established through observations. 
Extragalactic clumps of gas with mass densities of ~ 3 x 10"^"^ g cm~^ and a values in 
the range of 1 to 10 (Mq yr^^Y^'^, can probably not be detected with currently available 
instrumentation, but should be searched for, once those capabilities become available. 

Values of a ;^ 1 to 10 in the indicated units also make sense when we recall that 
a = e~^/^, as discussed in section 4. As a coefficient for star formation efficiency 
eMj , one can only expec t e to exhibit va l ues e < 1. If a is as high as 10, e would 
have to be of order 0.03. iKrumholz et al.l (120121) estimate star formation rates of 
SFR = /(if2)e//[32GpV37r]V^ where /(i/a) is the fraction of gas in the form of molecular 
hydrogen. They estimate the free-fall efficiency /(if2)e// to be of order 1%. Within 
reasonable uncertainties, this suggests that an estimate of e ~ 0.03, and thus a ~ 10 also 
lies within an acceptable range. 
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8. Summciry and Discussion 

We have shown that for a large fraction of the galaxies whose central regions were 
observed in the Sloan Digital Sky Survey, the observed gas metaUicities, Z^., are consistent 
with (i) a quantitative physical relation for star formation through episodic infall of gas 
devoid of heavy elements; (ii) thorough mixing of infalling and native gas before onset of 
star formation; (iii) a star formation rate proportional to the 3/2 power of the infalling 
mass rate, Mj; and (iv) intermittent quiescent phases devoid of star formation during which 
the native gas in a galaxy exhibits a characteristic elevated gas metallicity, Zq, dependent 
on galaxy mass, M*, and a characteristic mean ratio of stellar mass to native mass of gas, 
Mg. Most if not all our SDSS star-forming galaxies with < 2.0 x IO^^Mq, and many 
with > 2.0 X IO^^Mq and large radii appear fed by infall. Sets of galaxies not obeying 
these relations invariably have masses > 2.0 x 10^°Mq and exhibit high metaUicities at 
high star formation rates, suggesting that for these galaxies star formation independent of 
mass infall plays a significant role. 

These features are largely reproduced by a mean-field-theory toy model for the physics 
of infall, which accounts for the (SFR) oc Mf^"^ relation in terms of the pressure infalling 
gas exerts on native gas in a galaxy to trigger, or at least accelerate, star formation. The 
model also permits us to estimate the mean densities and velocities of clumps of baryonic 
matter traversing the dark matter halos in which the SDSS galaxies may be embedded. 

This is as much conformance to our data as can be expected from a mean-field-theory 
toy model, which of necessity disregards variations in velocity V, variations in clump 
densities p, and selection effects that favor detection of larger, more massive galaxies at 
higher redshifts, and smaller, less massive and less luminous galaxies at lower redshifts. 
But it suggests that the basic features of the model capture the essence of population 
distributions and metaUicities found in SDSS star forming galaxies. 
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A number of papers publis hed in recent 



chemical enrichment in galaxies (IBrooks et al 



years 



lave discussed t 



2007 



le absolute value of 



Kobayashi et al. 



2007 



Finlator fc Dave 



20081 ) : modeling the processes involved requires consideration of the chemical constituents 
of supernova ejecta, their ejection velocities, the gravitational potentials of the parent 
galaxies, and thus calculations on the rate of retention of metals in galaxies of different 



types and masses. As pointed out by iGenel et al. 



fl2009h 



past merger history can a lso be 



Leitner fc Kravtsovl ( 120111 ) have 



an important factor. A number of authors, most recently 
also considered the extent to which winds from evolved stars may generate sufficient gas 
to keep up star formation within quiescent galaxies at a low level on a continual basis, 
particularly in massive galaxies. While of great interest, a detailed understanding of these 
effects will require calculations and detailed modelin g be yond the scope of the present 



paper, as do recent attempts by 



Bouche et al. 



mm and 



Shi et al. 



(1201 ih to model SFRs 



through new star-formation laws. Much of this theoretical work is still tentative and will 
need to be tested by accumulating observational evidence. Some theoretical predictions 
could potentially be tested against the quantitative metallicity relations brought out in the 
present paper. 

Of particular interest to such studies may be the upper limits to the metallicity of 
infalling material that our study indicates, ~ 0.125 x 10"^ ±0.125 x lO^^ N(0)/N(H). This 
appears to be a useful constraint on the metallicity of material transiting galaxy halos and 
falling into the galaxies. 
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